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ANALYSIS STABILITY OF THE NONLINEAR EPIDEMIC MODEL 
WITH TEMPORARY IMMUNITY 


LAID CHAHRAZED 


ABSTRACT. The present paper present a nonlinear mathematical model, which analyzes 
the spread and stability of the model epidemic. A population of size N(t) at time t, is 
divided into three sub classes, with N(t) = S(t) + I(t) + Q(t); where S(t),I(t), and Q(t) 
denote the sizes of the population susceptible to disease, and infectious members, quarantine 
members with the possibility of infection through temporary immunity, respectively. This 
paper deals with the equilibrium and stability, precisely the global asymptotic stability of 
endemic equilibrium under certain conditions on the parameter. Second stochastic stability 


and finally the equilibrium and stability of the epidemic model with age. 


1. INTRODUCTION 


Motivated by study the dynamics of the different nonlinear epidemic models in all the 
references, the present paper addresses a mathematical model, which examine the evolution 
of the infection, reported in the total population. 

The spread of the epidemic might naturally dependent on possible contacts between sus- 
ceptible and infectious. 

In this paper, we discuss the equilibrium and stability of the model with temporary im- 
munity and the different positives parameters. We have made the following contributions: 


e The equilibrium and stability of the model, we obtain a disease-free equilibrium in 
the absence of infection but in the presence of infection, it was a unique positive 
endemic equilibrium and we define the basic reproduction number of the infection . 

e Second the global asymptotic stability of endemic equilibrium. 

e Next, we introduce a Brownian motion to system and we transform it into an Itó 
stochastic differential equation. 

e Finally study the equilibrium of epidemic model with age. 


Key words and phrases. Basic reproduction number; endemic equilibrium; global asymptotic stability; 


epidemic model; lyapunov functional; temporary immunity; stochastic stability. 
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2. MATHEMATICAL MODEL 


This paper considers the following epidemic model with temporary immunity: 


SHS : + pp — vy — (m +d) S(t) — BS(t)Q(), 
(1) Í (t) = BS(t)Q(t) — (us + d) I(t) — ye S(t — 7)Q(t — 7), 
Q (t) = ye?" S(t — 7) Q(t — 7) — (us + d) Q(t), 


Consider a population of size N(t) at time t, this population is divided into three subclasses, 
with N(t) = S(t) 4- I(t) + Q(t); where S(t), I(t), and Q(t) denote the sizes of the population 
susceptible to disease, and infectious members, quarantine members with the possibility of 





infection through temporary immunity, respectively. The positive constants 44, Hə, and 
l3 represent the death rates of susceptible, infectious and quarantine. Biologically, it is 
natural to assume that u, < min {4,43}. The positive constant d is natural mortality 
rate. The positive constants u represent rate of insidence. The positive constant y represent 
the recovery rate of infection. The positive constant 8 is the average numbers of contacts 
infective for S and J. p the positive constant is the parameter of immigration.v the positive 
constant is the parameter of emigration. The term ye "?* S(t — T)Q(t — T) reflects the fact 
that an individual has recovered from infection and still are alive after infectious period 7, 
where 7 is the length of immunity period. 

The initial condition of (1) is given as. 


(2) S (n) = 9i(p, 1) = 92(, Q (n) =n), | —- TS SO, 


Where ® = (1, ®2, ®3)T € C such that S (y) = 4$4(y) = ®,(0) > 0, (n) = ®2(n) = 
$5 (0) = 0, Q (n) = $5 (n) = $3 (0) = 0. 

Let C denote the Banach space C([—7, 0], IR?) of continuous functions mapping the interval 
[77,0] into IR?. With a biological meaning, we further assume that ®; (n) = 9;(0) > 0 for 
SESS 

Consider the system without the parameter of emigrations. Hence system (1) can be 


























rewritten as 
S (t) T i (t) - BS (t) Q (t), 
(3) I (t) = 8S (t) Q (t) — (us + d) I (t) - ye ^" S(t — 7)Q(t — 7), 
Q (t) = ye S(t — 7)Q(t — 7) — (ug + d) Q (2), 
With the initial conditions. 


(4) S(n) = ®1(n), (n) = 92(n) ,Q(n) = ®3 (n), =F <N <0, 


Where 6, (0) > 0, 65 (0) > 0, 65 (0) > 0, —T<7~0. 
The region Q = {(S,J,Q) € R$, S+I+Q<N< ol is positively invariant set of 
(3). 
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3. EQUILIBRIUM AND STABILITY 


An equilibrium point of system (3) satisfies 
ptpu-—v—(u-d)5 —85Q — 0, 
(5) BkSQ — (Ha + d) I — ye "^7 S(t — r)Q(t — 7) — 0, 
ye He S(t — T)Q(t — T) — (ug + d) Q — 0, 
We calculate the points of equilibrium in the absence and presence of infection. 
In the absence of infection J = 0, the system (5) has a disease-free equilibrium Ep: 


Lo ae eiui = » 
y= (8 1, 6) = (SFE, 0.0). 
p, +d 


The eigenvalues can be determined by solving the characteristic equation of the linearization 
of (3) near Ep is 





—(f,+d)—A 0 — Sele 
(6) et 0 -Qs*d-A SANG |-o 
0 0 ee c (ud) cd 


So the eigenvalues are 
Ai — — (ui t d), Ao — — (uo t d). 


In order for A1, A», to be negative, it is required that. 





Jd v)ye "n 
(7) (p+ u — v) 
p, +d 


Then we define the basic reproduction number of the infection Ro as follows. 


< ua +d 


(0 u ye ar 
(u + d) (us + d) 


In the presence of infection J Z 0, substituting in the system, Q also contains a unique 


(8) R= 





positive, endemic equilibrium E* = (S*, I*, Q*)’ where 


ko ptu—-v d. 
St (0 pad "s Ro? 














(9) Dd) [edger — niin 
Qt = SF (Ro — 1) 
Note that 
x (pt p—-v)(u t) p+u-Y 

NX = | 10 
Cae eee ~ 

| (Ro = 1) (m + d) (Lo = lg) 

p (Lo + d) 


So E* = (S*, IT", Q*)" is the unique positive endemic equilibrium point which exists if 
Ro > 1. 
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Theorem 1. The disease-free equilibrium Eo is locally asymptotically stable if Ro < 1 and 
unstable if Ro > 1. 


Theorem 2. With Ro > 1, system (3) has a unique non-trivial equilibrium E*is locally 
asymptotically stable. 


4. GLOBAL ASYMPTOTIC STABILITY OF ENDEMIC EQUILIBRIUM 


Consider system (3), with introducing the variables, 


z(t) - S(t) - S5, — wt)—I(t) - I, — z(t) - Q(t) - Q7, 
System (3) is centered at the endemic equilibrium Et = (Sz, I*, 02)”, then 


a(t) = [- (44 + d) - 8Qz] z  [- 881] z, 
(11) g(t) = (8 — ye) Qt] [7 Qua + d] v + (8 — 7") S$] z, 
a(t) = [bye Qt] x + [ye 82 — (us + d)] 2 


Lemma 3. Let 
S (s) 2 S(0) >0, I (s) 2 I (0) > 0 for all s € [77,0] and Q (0) > 0. 
S(t), I(t) and Q(t) solutions of system (3) are positive for all t > 0. 


Proof. For contraduction there exists the first time to, such that S (to) Q (to) = 0. 


e Assume that S (to) = 0, then Q (t) > 0 for all t € [0, t9]. With Eq 1 in the system (1) 
we have 
S (tg) - p - un —v » 0. 
For S (to) = 0, So > 0, we must have S (to) < 0 which is contradiction. 
e Assume that J (to) = 0, then with Eq 2 in the system (1) we have 


I (to) = —ye ^" S(t — 7)Q(t — 7) 


I (to) is positive because S(t) and Q(t) solutions of system (1) 
are positive for all t > 0. 
e For J (to) = 0, J > 0, we must have Í (to) < 0 which is contradiction. 
e Assume that Q (to) = 0, then S (t) > 0 for all t € [0, to]. with Eq 3 in the system (3) 
we have 


Ò (to) = ye™ S(t- 7)Q( — 7), 
Q (to) = y f e«-9s(s)Q(s)as. 


to—T 


to 
S (s) > 0, S (s) > 0 for all t € [0,to). We have y f e-s»-9(s)Q(s)ds > 0, and 
to—T 


Q (to) — 0, which is contradiction. 
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Lemma 4. Let 
S (s) = So > 0, Q (s) = Qo > 0 for all s € [—7,0] and Qo > 0. 


Then 
pHH=VY 
p, +d 





50 < max f rr Qo} =M 


Proof. We have 
N(t) = S(t) + I(t) + Q(t), 
For Ro < 1 the solutions S(t), I(t) and Q(t) approach the disease free equilibrium as 


t — oo. 
With Eq 2 in the system (3) we have I € — (u + d) I, hence if ji; +d < 0, 


lim I(t) = 0, 
t—oo 

With Eq 3 in the system (3) we have. 
lim Q(t) — 0, 
too 


With Eq 1 in the system (3) we obtain $ = p + u — v — (p, + d) S. 


lim S(t) = BERE SE 

t>00 I +d 
Hence. 

lim N(t) = AREA 

too My + d 


From lemmal, S(t), I(t) and Q(t) solutions of system (1) are positive. 


S(t) < ENE m allt < 0. 
p, +d 


Suppose that 
ptpu-v ptu-v 
N(0) < serra then N(t) € ————— 


Hy + — py +d 
On the contrary 


If N(0) > S STA then N(t) < N(0), and S(t) < N(0) for all t > 0. 














Theorem 5. Let S(s) = So > 0, Q(s) = Qo > 0 for all s € [—7,0] and Qo > 0. EX is 
globally asymptotically stable for all T 








1 In wM+3wQ* 
y 2w(n, +d) ? 
1] wM+3wQ* 
n 
T > MAI Xs dus rd)-BM | 
1 In Qe —3wM 





3 ^ X(ug*d)-BM 
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Where 
prey 
py +d 
8Q7 
My + pi + 2d 





M = mad 


Proof. We consider system (3). 
Let us introduce the functional 


V (53,2) = go (oo) e 5 (2), 





The derivative V (x, y, z) is 


V (x,y,z) = wit tyu) (THU) + yy + 22 
(— (u + d) — BQ*) x — BS*z + (B — ye "27) Q*x 
= (u td) y+ (8 — Ve") Siz 
+y [(8 — ye”) Q?x — (ma + d) y + (B — ye”) S22] + 
«s [Bre Qa + (16S — (s + d) 5 
= -u(u +d)? — |(w + 1) (us +a] y? 
—(ye e = (ia e 
+[8Q} — w (u + d) — w (My + d)] xy 
yana KA aa t=T) 
— (w +1) 9 ye — 7) 
GQ) ye t zet — T) wanye Pp) 


— (w + 1) Siye? yz(t — T) + Szye #7 zz(t — T). 








= w(x+y) 











By lemma 2 we have S(t) < M for all t > 0 and w is an arbitrary real constant choosing 
as follows 


TERES. HR 
Hi + Ha + 2d 
V (2,92) € —w(m-d)z^— [(w- 1) (m  d)] y’ 

— (us + d) 2? 
+BMyz — [wQ*ye"2"| zx(t — 7) 
- (w +1) Qie" yx(t — T) 
+Qžye 2 za(t — T) —wM ye?" zz(t — T) 
— (w + 1) Mye" yz(t — T) 
+Mye "7 zz(t — T). 
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Applying Cauchy-Chwartz inequality; we obtain: 


V(z,z) € —w(u+d)2?—-[Wt]1)(mM+d]¥ 
— (us + d) z 
1 
-3uQtye n |a? + a? (t — 7)] 
1 
-1 (w +1) Qe fy?  z'( — 7) 
1 * — HoT 
*3Qre HaT [2? + z^ (t — 7)] 
~ SM ye" [z? + 2?(¢ — 7)] 
1 
E (w +1) Mye 2? b? + 2?(t — 7)] 
1 
*t9Mwe [2? + 2*(t — 7)| 
1 
18M s 2], 





IA 


[aie 4) - genet (M+ Q2) | 2 


5M - (w +1) (ja +a) - 56 Doe (M+ Q2)] v 





BM ~ (p +d) + je (M e | 2 


- leQtye ^7] z^(t- 7) — E (3w + 1) Met" z2 (t-r) 


Choose the Lyapunov functional 


t 
V (24, Ut; 24) = V (x, Y, z) i uu e rat J a^ (0)d8 
i—T 
1 t 
x (3w - 1) yMe^"7 f z?(8)d0 


i—T 
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Then 


V (T4, Yt, zt) = V (7, Y, z) = wQ*ye 27 x? (t) 
t+wQ*ye "Tr (t — 7) 
—- (3w + 1) yMe ^ z?(t) 





RE (3w + 1) yMe 7 z?(t — 7) 


IA 


7 (n +d) — zoe (M + 22) E 


Eos (w + 1) (uS + d) 
i (w +1) ye" (M + Q1) 


2 








T [5^ BM — (ug + d) + ec (M + Q% |? 
— leQtvye "7] z? + wanye "72? (t — 7) 


1 
Sh 3w + 1)yMe7 zi za 
1 


+5 (3w + 1) yMe "?* z?(t — 7), 


IA 


[aie es d) - gene m (a Qi 


+ [58M - (0o - 5o tre ar Q2] 0? 





+ EU (ua +d) + see (M +Q% | 


1 
- [wQtye "7] z? — É (3w + 1) ymer 2 
Therefore 


Vistas ied l^ (i, d) Sere" (M + 97) z? 
(o +1) EE - 48M 
+z (w + 1) ye 7 (M + Q7) 


= (us + d) — 38M E 
-3ye "(Qo — 30M) |^ 


y? 
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While the above inequality is always negative provided that 


wM+3uQt 
2w (p14 +d) 


1 l wM+3wQ7 
n 
T 2 max4 a D e ora oM 
1 ln Qž—3wM 
y T 3us+d)—8M 


Ln 
- 








With application of the Lyapunov-LaSalle type theorem in [10] 


lim z(£) 20, limy(t)=0, lim z(t) — 0. 


ioo i—oo ioo 














5. STOCHASTIC STABILITY 


We limit ourselves here to perturbing only the contact rate so we replace 3 by 8 -- oW (t), 
whereW (t) is white noise (Brownian motion). The system (3) is transformed to the fol- 
lowing Ito stochastic differential equations, with ^y = »ye "27 


dS =|[p+p—v— (m + d) S — BSQ] — oSQdW, 
(12) dI = [BSQ — (us + d) I — y95Q] + cSQdW, 
dQ = [y98Q — (us + d) Q], 
In this section, we will proof, under some conditions, that Eo is globally exponentially 


mean square and almost surely stable, and for this purpose, we need the following Theorem 


Theorem 6. The set Q is almost surely invariant by the stochastic system (12). Thus if 
(So, Io, Qo) EQ) j then P (S, 1, Q) € Q| zu 


Proof. The system (12) implies that dN < [p + u — v — (u + d) N] dt, then we have 


ptu-v p 
N (t) € ————— No — ——— |, for all t > 0. 
(t) < m «(s ze) or allt > 


Since (So, Io, Qo) e Q, then 
uci tu-—v 
— mtd 


There exist €9 > 0, such that So > £o > 0, Ip > £o > 0 and Qo > £o > 0. 
Considering 


(13) N (t) , for all t > 0. 


ve = inf{t>0, S(t) «eor I(t) € eor Q(t) € e£), fore € £o, (14) 
Ors lim ve = inf {t > 0, S(t) < 0 or T (t) < 0 or Q (t) <0,$ 

Let 

ptu-v 

(u + d) Q (t) 





V (t) = log EE end EUN uid 


(iso Cura S 
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Then, using Itó formula we have, for all t > 0 and T € [0,t ^ vel, 

















_ ptu-v 
dV(T) = dT 4-cQdW 
mw m 
(vo — 8) F SQ 
+ dT + o —dW 
+ (u + d) + 18? DE 
t [7*98 + (us + d)] aT, 
Hı + Ha + H3 + 3d Q 
< zw p IL 
(15) dV (T) < -8Q - 1I 1s? dT +07 (I — S)daW 














With (13), we have S,J and Q € D uu 








(4 +d) 
Let 
L = py + My + us + 3d (16) 
peta | II 
| p, +d p, +d 
Q 
I = > 
fü) P 


We remplace (16) into (15) , we obtain 


(17) dV (T) € LdT + o (I (T) — S (TY) f (1 (TY) aW, 
Then 
(18) V(T) & LT 4 ef (1 (u)) E (u) — 8 (u)) dW (u), 


p 
With proposition 7.6 in [6], o f f (1 (u)) (1 (u) — S (u)) aW (u) is mean zero process then, 
0 


(19) EQ Cypr 


for allt > 0 and T € [0,tA vel, 





S(tAv),I(t ^v), and Q(t^v) € E P i 


(m + d) 
Then 
E(Q (t ^ve)) < BAV) < Lt, 
V (tAv) = 0, 
(20) E(V (^v) 2 E(V (t)) x huey 2 P (ve < t)log PET 


Where kj,.<4 is the indicator function of a subset [ve < t], 
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Combining (19) , and (20) , we obtain 


(21) P(w-<t) <a for allt>0; 





Drs 


for all t > 0, and € > 0, we obtain P (v < t) = 0; 
From where 





Puro) e 














6. THE MODEL WITH AGE 


The age distributions of the numbers in the classes are denoted by S(a,t), I(a,t), and 
Q(a,t), denote the sizes of the population susceptible to disease, and infectious members, 
quarantine members with the possibility of infection through temporary immunity, respec- 
tively of age a, at time t, d(a) is the age-specific death rate, 

The system of partial equations for the age distributions is 

ot Se = — (us + d(a)) S(a, t) + 8,(t)5(a, t), 
(22) Gr $2 = Fi (t)S (ast) — (us + d(a)) (a t) + v (t — T)S (ast — 7), 
e e = -q(t - 7)S(a,t — 7) — (us  d(a)) Q(a, t), 
With 
Bit) = —BfQ(a.t)da (23) 
mtr) = —vfe™"Q(a,t — 7)da 


6.1. Equilibrium and stability. Assume that sub population does not depend on the time 
when the system (22) is written as follows 


4$ — (8, — i4, — d(a)) S(a), 
(24) al — (yi — 83) S(a) — (us + d(a)) I(a), 
49 — —YyS(a) — (us + d(a)) Q(a), 


The initial condition of (24) is given as 
(25) SO)=S:, 1(0-h, Q(0)=Q1 


Differential equations of the system (24) are solved with different methods of resolutions and 
with (25), so 





(26) S(a) = S e7 Hoa $ (a), 
(27) I(a) = I, ®(a)e "?* E ——À M (e- 67600 = gat) 
(28) Q(a) = Qi P(a)e “3% uM Sla) | (e (7609 E ena) 


Hı — By — ps 
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Where 
(29) (a) = exp (— f d(a)da) 
The system (24) has the unique positive equilibrium point P}, 
Roo oU z 
P, = (8, I 01) = (0, 0, 0)”. 


We calculate the Jacobian matrix according to the system (24) with Pi 


bı — u, — d(a) 0 0 
J (Pi) = À — Y — (us + d(a)) 0 
—Yo 0 — (us + d(a)) 


The epidemic is locally asymptotically stable if and only if all eigenvalues of the Jacobian 
matrix J(P,) have negative real part. The eigenvalues can be determined by solving the 
characteristic equation of the linearization of (25) near Pj is 


bı — m — dla) — A 0 0 
(31) det À — Yo — (u + d(a)) — A 0 —0 
= 0 — (us + d(a)) - A 


So the eigenvalues are 
Ai = fı =p —d(a), A= — (us + d(a)), As = — (us + d(a)) 
In order to A, As,and As; will be negative, it is required that 
By < us + d(a) 


The basic reproduction number A is defined as the total number of infected popula- 
tion in the resulting sub-infected population where almost all of the uninfected. The basic 
reproduction number of the infection Ro is defined as follows: 





By 
32 Ry = ——— 
Pd = uda) 
The time during which people remain infective is defined as 
1 
T——————— 
hı + d(a) 
The doubling time tg of the epidemic can be obtained as 
(In 2) T 
33 ta = 
(33) uer 


Theorem 1 The disease-free equilibrium P) is locally asymptotically stable if Ro < 1 and 
unstable if Ro 1. 
Let (26), so if Ro < 1 then u — 8, > 0, so S(a) converges to zero. 
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Let (27) , so 


(34) 


(51 — b1) $19 (a) 
Hy — By — be 





I(a) = now = | lace m; = min {Hy E Bi, Ha} 


If Ro < 1, i(a) converges to zero. 
Let (28) , so 


(35) 


71151P(a) 


Q(a) = [aiao - 2520 


| ems, m, = min {H, — b1, H3} 


If Ro < 1, Q(a) converges to zero. 


Conclusion 7. This paper addresses a epidemic model with temporary immunity, whenever 


the 


quarantine individuals will return to the susceptible. The endemic equilibrium is globally 


asymptotically stable, then under some conditions, study the stochastic stability. Finally, the 


equilibrium and stability of the epidemic model with age. 
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